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Abstract 

Starting from the von Neumann equation, we construct the quantum evolution 
equation for the effective action for systems in mixed states. This allows us to find 
the hierarchy of equations which describe the time evolution of equal time correlators. 
The method is applied explicitly to a scalar theory with quartic self-interaction. 
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Introduction 



During the last years, great progress has been made in the study of statisti- 
cal systems of quantized fields in thermal equilibrium. Many physically interesting 
quantities, however, are typically related to out-of equilibrium scenarios, like density 
fluctuations in the early universe, the radiation signature of hot quark-gluon matter 
or baryon number violating processes. The study of such quantities in systems which 
are not near to equilibrium requires a framework which is able to deal with mixed 
states characterized by general density matrices. 

Although the Schrodinger equation already determines the exact evolution of 
initial density matrices and can be integrated to solve the problem of time evolution, 
it is of less use for practical purposes. What one is interested in and what is used 
to describe a physical system is not the density matrix but a number of observables 
derived from it. At first sight it may appear to be rather a matter of taste to 
study time evolution of observables instead of the density matrix, but there is also a 
profound physical reason to change the point of view. Time integration of the density 
matrix requires the exact knowledge of the initial state, which is not accessible since 
observation is usually restricted to only a small number data. Consequently, the 
description of a system requires apart from the knowledge of the exact motion of 
observables also a statistical assumption replacing the incomplete knowledge of the 
initial state. 

In this letter we focus on the first step of the problem. We investigate on the 
exact counterpart of what the Schrodinger equation is for the density matrix: A 
quantum Liouville equation for the effective action which encodes in a compact form 
the time evolution of all equal time correlators being one basis of a complete set of 
observables. 

Construction 

The following construction will be formulated in the Schrodinger picture, that 
means operators are time independent and time evolution is contained in the density 
matrix. The resulting evolution equation, however, is independent of the picture 
chosen. 

Starting point is the von Neumann equation for the density matrix operator, 

id t p(t) = [H,p(t)}, (1) 

with Hamiltonian H. It is a functional of the position operator 3>(x) and its canonical 
momentum operator n(x). In a given mixed state described by a density matrix p(t), 
operator expectation values are given by the trace average 

(0) t = tv(Op(t)), (2) 

where the hermitian density matrix is normalized to unity, trp(i) = 1, which is 
consistent with unitary time evolution. 
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We introduce the generating functional of polynomials in momentum and position 
operators by 

Z[j, k; t] = (D[j, k]) t = tv(D[j, k]p(t)) (3) 

with current operator 

D\j, k] = exp(ij ■ $ + ik ■ II) (4) 

and shorthand notation J <i 3 xj(x) < I > (x). 

Since neither position nor momentum operators explicitly depend on time, inte- 
gration is done with respect to space only, and also there is no need to introduce time 
dependent currents. This is a very special feature of the Schrodiger picture where all 
the dynamics is hidden in the states and the density matrix. The definition breaks 
Lorentz invariance. However, it is broken anyway since in the Hamilton formulation 
the time direction (evolution direction) orthogonal to the time slice of quantisation 
plays a privileged role. 

The current operator can be rewritten as 

i % 
D[j, k] = exp(ij • $) exp(zfc • n) exp(-j • k) = exp(ifc • n) exp(r/ • <£) exp(— —j ■ k) (5) 

using the fundamental commutator [3>(x), II(y)] = i<5 3 (x — y) and the first non-trivial 
term in the Campbell-Baker-Hausdorff formula. From that representation, powers 
of position and momentum operators can be obtained by the functional differential 
operators, 

J- - ifc) D\j, k) = <Z>D[j, k], (Jj- + ij) D\j, k] = UD[j, k], (6) 
or, equivalently, by acting from the left with 

D[j, k] ^ + \k^ = D[j, D[j, k] ^ - ijj = D[j, k]U. (7) 

Polynomial operator expectation values 0[<&, H] can thus be cast by firstly construct- 
ing the operator with inverse order by reading it from the right to the left, secondly 
the replacements — > (5/i5j — k/2) and n — > {5/i5k + j/2) in the inverse ordered 
operator, third acting with this expression on the generating functional, and finally 
letting j, k — > afterwards. We note that the functional differential operator replace- 
ments have the negative commutator with respect to their field operator counterparts. 
We also note that if we used the representations $ — ► 5/i5j and n — ► 5/iSk instead 
of those of Eq. (||), we get the symmetrically ordered operator polynomial. 

Let us construct the time evolution of the generating functional. From its defini- 
tion and the equation of motion for the density matrix, we get 

d t Z\j,k;t]=itv(H[n,^]D[j,k]p(t))-itr(D\j,k]H[U,^p(t)). (8) 

In view of the arguments given above, the Hamiltonian can now be replaced by 
a functional differential operator which can be drawn out of the trace. We get a 
Liouville equation for the generating functional 

d t Z\j,k;t]=L[j,k]Z[j,k;t] (9) 
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with Liouvillian 

L fa -, *, = iH^± + ij, ^ - - ur^ - i,, A + i*]. (io) 

It is interesting to restore powers of ft at this stage. We note that the terms with 
factors one-half come from the fundamental commutator and are of relative order fi. 
Expanding the r.h.s. in h, the leading term cancels due to the relative minus sign in 
the summands in ([U]). The first non- vanishing contribution is linear in fi as the time 
derivative which also contains one h. That contributions correspond to the classical 
limit of the Liouville equation. 

Z[j, k; t] is the generating functional of all equal-time correlators corresponding to 
Green functions at equal time. The generating functional of the connected correlators 
W\j, k; t], however, is related to Z\j, k; t] by 

e wVM = Z \j, k;t}. (11) 

This relation implies for the r.h.s. of the Liouville equation 

L[j, k]e w[j ' k ' t] = L[W] [j, k; t]e w[j ' k ' t] (12) 

which defines a non-linear evolution operator. The exponential also appears in the 
time derivative of the Liouville equation and thus cancels out. Evolution is now given 
by the Liouvillian-like equation 

d t W[j,k;t]=L[W}[j,k;t] (13) 

This equation qualitatively differs from the counterpart for Z[j, k; t] since the operator 
L is non-linear in its action on W. 

In order to establish a relation to observables, we introduce the classical c-number 
field and momentum variables. We define the quantities 

5W[j,k;t] SW[j,k;t] 

which coincide with the expectation values ($(x)) i and (II(x)) i when the currents 
vanish. We want to point out that the classical variables ip(t, XL)\j=k=o and 7r(t, x)\j = k=o 
do depend on time and are independent of the picture chosen - in the opposite to 
their Schrodinger operator counterparts. In the definition (|l4|), the classical field and 
momentum are functionals of the currents j, k. Suppose we can invert the relation 
and express the currents as functionals j = j[tp,Tr], k = k[cp, tt], then the effective 
action defined by 

T[ip,Tr;t} = <p(t)-j + n(t)-k + iW[j,k;t] (15) 

can be expressed as a functional of the classical variables. It obeys a relation dual to 
the definition (|l|) 

§r[<p(t)MtY,t] _ . <STMt),7r(t);t] _ 
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That implies that the currents vanish at the minimum of the effective action where 
the classical quantities have their interpretation as operator expectation values. 

In order to apply terms correctly, let us remind the reader that the quantity 
r discussed here differs slightly from the definition used in standard field theory, 
namely, it corresponds to the generating functional of vertex functions at equal times. 
Here, the effective action itself has explicit time dependence from the classical fields, 
and an implicit one inferred from the evolution of the density matrix. The total time 
derivative thus behaves as dT/dt = tp-5T /5(p + Tr -5T /5ir + dfT, which, when compared 
with the time derivative of the definition (|l5l), implies 

d t Y\w,t]=id t W\j,k;t] (17) 

where we took into account the time independence of the currents and the relation 
(16). It is now straight-forward to translate the equation of motion for the generating 
functional of the connected Green functions into an equation for V. One only has to 
express the Liouvillian in the evolution equation ( |l3| ) in terms of tp and tt. In X[W], 
currents can be replaced by the relations (|l6|). The first derivatives of W are given 
in Eqs. (0), and higher derivatives can be found by successive differentiation of that 
relation. We get for the first few terms 

- iW Jj = <P,j = - tVV 

= — r jW (r t(pvv ,r )tp(p + r i(p(p7r r^ n 



-iWjkk = -r~i (r W7r r~^ + r ^r^) (is) 

where we wrote shorthand 5T/5ip =: T >(p etc., and W itself has to be replaced by 
— i(T — ip ■ T & — 7r • r jVr ). Along that line of arguments, we can translate the operator 
L[W] and finally get the desired Liouville equation for the effective action, to wit 

dt T = F[T]=iL[W]\ jMMipi7T] . (19) 

The classical limit of it was studied in |l| using a different approach. The relation 



(1£) is local in the time parameter and the determining equation for T. Expanding 



the effective action in powers of ip and tt 



n=l ' \j=l ' 



■ ••x n ;iMxi,i) • ..<p(xj,t)n(yL j+1 ,t) . ..7r(x n ,t) (20) 

3=0 



and comparing the corresponding coefficients at both sides of (19), one gets an infi- 



nite non-linear and highly coupled hierarchy for the time dependence of the proper 
vertex functions T^. The solution involves an infinite number of initial conditions 
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corresponding to the infinite set of observables necessary to describe a system com- 
pletely. Suppose we can find a solution of that hierarchy, at least a perturbative 
one, one may ask for the time evolution of observables. We know that the c-number 
fields ip and ir have their interpretation as expectation values of the position and 
momentum operator when the currents vanish. Thus, from the dual relation ( |i~6|) we 
find the effective equations of motion 

ar[y(*),7r(t) ; t] = 6r[<p(t)Mt);t] = 

6tp(t) ' 5ir(t) ' 1 ' 

as variational minimum of the effective action. These equations determine the time 
evolution of the average field (3>) t and momentum (II) t resp. Expectation values of 
symmetrically ordered higher order operators are given by the derivatives of W in 
Eq. fll8| ) evaluated at the solution of (pip. 



Scalar field 

Let us consider the non-interacting motion of a scalar field characterized by the 
Hamiltonian 

H° = ^J d 3 x (n 2 + $(-A + m 2 )$) (22) 
which, by Eq. (001), gives rise to the free Liouvillian 



Since, in this simple case, the operator is even linear in the functional derivatives, 
its action on the exponential is also linear in W, L°e w = (L°W)e w . Rewriting 
everything in terms of the effective action yields 

ar = F o [r] = j A ^ (X)( _ A + m2) JL_ _ , (x) _iL_) m 

which has a remarkably homogeneity property. The r.h.s. is a map of the class of 
polynomials in ip and ir of a given fixed order onto itself and the n-point vertices T* n 



in Eq. (20) obey a closed equation each. Time evolution can best be made explicit 



by a diagonalisation procedure using the normal coordinates 



^ = \/^ + ^= 7r ' w k = >/m 2 + k 2 , (25) 
V ^ \J Lio 



where the expansion coefficients in ( |20[) transform accordingly. The Liouville operator 
now takes the form 

^ [rl = i /A„( ? M_ r M) (26) 
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which implies for time evolution (the index j now refers to the number of £'s in the 
expansion in normal coordinates) 



k n ;0)exp(i^o; ki t - i ^ ^i). ( 27 ) 



=3+1 



As expected, each proper mode just oscillates harmonically. We want to point out 
that motion is perfectly compatible with arbitrarily high order vertex functions, and 
there is no physical argument to neglect them in the first place. In that sense, in a 
non-interacting system higher correlation functions have to be taken into account in 
order to allow for general initial conditions. 

In contrast to the simple form of the time evolution of the effective action, the 
solution of the on-shell conditions (^) involves to resolve a non-local non-linear 
equation for the classical fields (p and tt even in the non-interacting theory. 

The inclusion of a quartic interaction term H = H° + XJ d 3 x(p 4 modifies the 
Liouvillian by 



L (A) 



iX d 6 x\A 



fc(x) + 



i8j(x 



■k 3 (x) 



(28) 



This gives rise to the non-linear Liouvillian 



F X [T] 



X dx \ 4 



5T 



d 3 u J d 3 v J dW w) (x,u)G (w) (x, v)x 

5 3 r 



-3^(x)G (w) (x,x) + 9? 3 (x " 



+ G^)(x,w) 



5(p(u)5ip(v)5Tr(w) 




(29) 



The graphical representation of the complete time evolution including the free part 
is depicted in Fig. 1. The filled blob denotes the generating functional T, and full 
and dashed legs attached to it correspond to derivatives with respect to ip and tt 
respectively. The unfilled circle with one leg stands for a factor ip and with two legs 
for the two-point correlators G^^ and . The diamond stands for the free kernel 
operator —A + m 2 and the dashed circle for a factor ir. 

The first two graphs correspond to the classical non-interacting and the third 
through sixth graph to the classical interacting Liouvillian. Only the last graph is 
a quantum correction being of relative order fi 2 . The system can be expanded in 
powers of (p and tt which leads to a hierarchy of coupled equations for the vertex 
functions T J n . 

One may now start to attempt to solve the hierarchy for the vertex functions by 
cutting it at a certain order. However, this somehow pragmatic approach, however, is 
not guaranteed to yield a physically sensible result. Firstly, one runs into consistency 
problems with the hierarchy. Even worse, a check would require to calculate those 
terms which have been dropped in the cutting procedure in order to justify them 
to be neglected. Moreover, guided by the exact solution for the free case, where we 
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have already seen that higher correlators have to be treated at the same footing as all 
others, it is at first sight not evident why just the more complex dynamical structure 
should admit any such simplification. Secondly, cutting the hierarchy implicitly im- 
plies an assumption about the initial conditions for the higher vertex functions. It is 
not easy to see what physical initial preparation cutting the hierarchy and assuming 
particular initial values just for the lowest vertex functions represent. It is even more 
difficult to figure out what statistical assumptions about the system are hidden in a 
particular cutting procedure. 

To summarize, the Liouville equation derived here makes a statement about the 
exact evolution of a system in mixed states. Any statement about the thermodynamic 
properties of that system, and in particular about how it approaches an equilibrium 
state cannot be read of from the dynamics but requires the application of the sophis- 
ticated tools of statistical mechanics. This goes beyond the scope of this letter and 
will be presented elsewhere. 
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